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ABSTRACT 


The purpose of this analysis is to determine the effect of transient heat flux on heat 
flux sensor response and calibration. 

We have begun this work by studying a. one dimensional case. This was done to 
elucidate the key parameters and trends for the problem. It has the added advantage that 
the solutions to the governing equations can he obtained by analytical means. The 
analytical results obtained to date indicate that the transient response ol a heat flux sensor 
depends on the thermal boundary conditions, the geometry and the thermal properties of 
the sensor. In particular it was shown that, if the thermal diffusivily of the sensor is small, 
say for a material such as Nicliromc, then the transient behavior must be taken into 
account. 
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INTRODUCTION 


I lie purpose ol this analysis was to determine the effect of transient heat flux on 
sensor response and calibration. 

Consider the case ol a heat flux sensor mounted in a wall subjected to both radiative 
and convective heat, transfer on the top side and purely convective heat transfer on the 
bottom side as shown in Fig. 1, where x and r are the axial and radial coordinates, 
respectively; L and f, I) and d are the thicknesses and diameters of the two layers; k is the 
thermal conductivity, a is the diffusivity, T is temperature. Ip and h 2 are the convective 
heal, transfer coefficients, q is the radiant incident heat flux, and T is the ambient 
lemperature. The subscripts i and 2 refer to the different materials of the sensor. 


In this analysis we make the following simplifying assumptions 

I II neglect re— radiation; 

2) constant thermal properties; 

31 constant, uniform ambient temperatures; 

4) constant, uniform convection coefficients; 

51 constant, uniform radiant incident heat flux; 

6) perfect thermal contact between the regions of different k; 

(7) the equations arc one-diniensional; 

Ste ady - State Solution To determine a base line case we first solve Tor the 
steady-state case. 

For this case, the governing equations and corresponding boundary conditions are, 
under the above assumptions. 


( 2 ) 


(V 


d 2 T, _ n 
fix 2- ~ 0 

d 2 T 2 _ n 


0 < x > f 
t < x > L 


dT 


+ 1 ' 1 (T , .-T fl9l )-q 0 = 0 at x = 0 
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, dT, 
k| TTx 

, dTo 

= k - tlx J 

II 

(I) 

i dTo 

k2 fr.'T” 

_ = h, (T, -T o 

) at v x = L 

(• r >) 


(3-a) 

(5-a.) 

% 

where lj = T| + i 

oo N | 

r» = To 

“x 

The solutions to equs. ( I ) and (2) are 

Tj = At x + Bj 
To = A o x + Bo 

The coefficients A|, B|, Ao, Bo can he determined iti terms of the Iwmndary conditions 
(3-a). (1), and (o-a). 

't he sensor output is based on the temperature difference bet ween two fixed points 
where the two alumel wires are attached. Assume that the coordinates of these two points 
are xq and xo respectively (they are usually at different layers), then 
T, (x,) = A,x, + B, 

To (xo) = AoXo + Bo 

Hence, 

AT S = T, (x.)-To(xo) (fi) 


(d) and (5) are rearranged as: 

-k, I'i I'i - h,l, at x - 0 

ko — + hoTo = hofo at x = l 
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where 


A| — A_. , Hi — (k| A i + li| fj )/h i 

— B, + f (Ai — A 2 ) 

kikali i + I< i hj I 12 L + + k|k2li2 

Note that AT S is independent, of time. 


Unsteady — State Solution 
I 'or unsteady state case, we have 

n < = lr- o< x <o>o 

a 2 0 - = lr- 0 < x < l, t > 0 


;Tr 

~jj^r - + ii|T| = h 1 1 1 at x = 0, t > 0 


T , = T , 


at x = f, t. > 0 


- Atx=f ’ ,>o 
~ — h hoTo — IhIN at x = Tj, t, > 0 


T, (x,0) = T, (x,0) = T w (|0) 

I he solutions to crpis. (7) and (8) with iionhoniogcncous boundary conditions (0) and (10) 
can be expressed in the form 

Tj (x,t) - ft (x,t) 4- (x) b + V’j (x) f> (II) 

i = 1,2 

where the (unctions y?j(x), i/> j(x), and ft (x,t) are the solutions of the following subproblems: 
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!•>> 43 = 0 


<lx^ “ 0 

wit li k| — f l'i vi = l*i F| 

at x = 0 

V 7 ! ~ V 7 ! 

at. x = f. 

1. fly>l _ Cl <£> 

k * ar~ - ! ' 2 ar- 

at x = ( 

ks + h 2 T 2 = o 

at x = L 

Therefore, 


V 7 ! ( x ) = f~'l x + Oi 
<£> (x) = Cox + Dj 


<l'T. 


d>) fI x 2 = 0 


d 2 02 

dx 2 =0 


with — kj ^-1 — 1- It |0| = 0 

at x = o 

0, = 02 

at x -- ( 

k 'aF- = k 4^- 

at x = ( 

(1 i f 'i 

1<2 (|x + lljXj = 1 1 2 1 3 

at x = b 

we have 


0i (*) = f-i x + F, 


0_> (x) = Ivj X + I’j 


lii this case, the sensor output will be dependent on time. 

To determine the 


output we calculate the temperature difference to be 
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■AT 11 (!,) = T| (x,t)-T, (x,,t) 

= 0\ (X|,t) — 0>( X_>,t) + [ <^|(X|) — ( x 2 ) ] T | 

+ [V ; l( x l) _ V*i( x 2)]f2 

Computational results sliovv that 

[ < r 7 l( x l) — V > 2( x - , )]f| + (0l(X|) — = AT S 

Therefore, 

AT u (l) = 0 { (x u i) - 0 2 (x 2 1) + AT S 

(c) 0,(x,t) and $>(x,t) are determined from the following homogeneous 
problem: 

d 2 o x m_ 

n '<)x- = dt 
d~ O2 d 0 2 
02 ex ' 2 = d\ 

do 1 

Subject to -ki Ox + h|0j = 0 at. x - q 


II 

at x = i 

00 x 00-2 


k 1 dx = k , c)x 

at x = ( 

m 2 


ko dx + 1 1 2 0> = 0 

at x — b 


and /A(x,0) = T Q - ^(xjl, - ^(xjlh 
let 0j(x,l) = X j ( x ) r(t) 

=> r(t) = c _/j n 1 
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win ‘io /in are eigenvalues which will lx* determined later. 


x „. <*' “ A ,» Shl 1 ?=r *> + B i» 


i 

_r 

f V 

A, 


X. JU (x)«A„ n Sin (^=x) + B 


2 II 


v^T 


on 


Cos {—^ — x) 

0 n 

Cos ( ■■ - x) 


Without loss of generality, wo choose A|ii = I, B , A () , and h >n can bo determined with 
the corresponding boundary conditions(and several pages of algebra). 

Finally, 

AT u (t) = S C n [X |(1 (x ) - X 2I) (x 2 )]e^n*' + AT S (13) 

n=l 


where C is determined by the initial conditions. The equation for the determinat ion of 


the eigenvalues /? is 

sin 7 + Hi 7 sin 7 
K cos 7 - KHi7sin 7 


—sin t] 


COS )j 


0 Il 2 »; cos (jij) + sin(yi/) cos(j//)-U-»//sin(jj/) 


— (()S // 
—sin ) 1 


=0 


where 7 = ^n f , tj — ^ 


b f. 


vCTT 


vTo 


1 - K t 

'■-In? 


K - K| 


FT 


0| 


n > 


11 


v/ 


Ivj 

lid. 


Some additional alegebra will lead to the following transeendtal equation lor the 


determination of the eigenvalues /Ju: 
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•■Jit— _ /3n) 


i - In f , , 
K j— J7- tfil 


i,-c 


AT 


/Ai)- 


(M) 


'lewril o equal ion (13) as: 


AT n (t) = S C e /J " u + A r I 

n = l 11 


-pwh 


rriS 


(15) 


where 


C „ = C„|X lp (x,)-X a „'(x)] 

Numerical example shows that pr « pr « Pr << , so we can just lake the 

first term to approximate the resull. In other words, we have: 

ui -pth . A r„S ’ : • 


hut 


AT 11 = c,e -/jrt 4- AT 
AT” = 0 


at, t = 0 


Therefore, ei — — A1 


r„S 


A'l 


,u 


Then, 


AT S “ 


Pr t 


(Ifi) 


where pt can be determined from Equation (1-1). The following observations can be made 
by examining F.qs. (14) and (Hi): 


(1) AT u (t) = AT s 


MS t -» x ; 


(2) The transient process depends on ( he boundary conditions ( 1 1 1 J i > ) , 
the geometry (Lb) and the properties (k|.kj,ri|.n.i) of the sensor. 

f 1 

(3) The transient term (e 1 1 ) may be very important if p t is small. 
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Typical Kxnmple 

To illustrate tl it: impact of these results consider a typical example. Assume that D 
n.S cm, f = 0.02 cm, L = 0.11 cm, choose h 2 /hi = 1.0. See Fig. 1. The material of the 
sublayer is alumcl whose properties are k> = >17 w/m-k, = 9.7 x 10 nr/s. If the 
materials of the foil are pure nickel and nichrome (811% Ni, 20% Cr), respectively. The 

result s are shown in Table I . 

Table l 

Foil Material o, (nr/ s) lh e ■ '(at t = 0.05 second) 

Nickel 2.3 x 10 -5 16.745 S.15 x 1 0~ 7 

(pure) 

Nichrome 3.1 x I0~ ^ 6.138 1.26 x 10 ' 

(80% Ni, 20% Cr) 

As can be seen from Table I, if we choose Ni— Cr alloy as the foil material, the 
—3 - 1 

transient term "e 1 " would be important. 

SUMMARY 

YVe have initialed a study of the transient heat flux response of heat I lux sensors. 

I his paper presents a one dimensional analysis of the problem. This was done to elucidate 
the key parameters and trend for the problem. The results of a heat llux sensor depends on 
l he thermal boundary conditions, the geometry and the thermal properties of the sensor. 

In particular, it was shown that if the thermal dilfusivity of the sensor is small, say for a 
material like Nichrome, then the transient behavior must be taken into account. 
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Schematic of Transient Analysis 



Fig. 1 
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